Convection in rotating annuli: Ginzburg-Landau equations with tunable coefficients 
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The coefficients of tfie complex Ginzburg-Landau equations tliat describe weakly nonlinear con- 
vection in a large rotating annulus are calculated for a range of Prandtl numbers a. For fluids with 
(J ~ 0.15, we show that the rotation rate can tune the coefficients of the corresponding amplitude 
equations from regimes where coherent patterns prevail to regimes of spatio-temporal chaos. 
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The complex Ginzburg-Landau equation (CGLE) 

dtA^A+{\ + ici)dlA - (1 - ic^)\A\''A , (1) 

which describes slow modulations of an ampUtude or 
envelope A near a Hopf bifurcation in spatially ex- 
tended systems, has been used extensively both to study 
nonequilibrium pattern formation flpl and as a model 
system for spatio-temporal chaos The qualitative 

dynamical behavior of solutions of the CGLE depends 
on the coefhcients ci and C3, which, for a given system, 
can be obtained from the underlying equations by labori- 
ous calculations; these have been carried out for a small 
number of systems only For ci and C3 small, as well 
as close to the line C3 — — ci, the dynamics is close to 
that found in the relaxational limit ci = C3 = 0, whereas 
for |ci| and |c3| large, the CGLE reduces to the nonhn- 
ear Schrodinger equation. In recent years, the compli- 
cated and often surprising dynamics that occurs away 
from these limits, has been intensively studied [^ji). In 
particular, it has become clear that the CGLE shows 
various regimes of spatio-temporal chaos when C1C3 > 1 
— see Fig. 1. The precise nature of the various chaotic 
regimes, as well as the existence and nature of the tran- 
sitions between them, is still under active investigation 
in the field of spatio-temporal chaos. 

In order to be able to investigate these chaotic regimes 
experimentally, one would like to have a system where 
the coefficients of the corresponding CGLE can be tuned 
in a convenient way through the spatio-temporal chaotic 
regimes. When it was discovered that a forward Hopf 
bifurcation to a quasi-one-dimensional wall-mode occurs 
in rotating Rayleigh-Benard convection in bounded con- 
tainers it was realized that the rotation rate might 
serve to tune the coefficients of the corresponding CGLE 
0. Kuo and Cross therefore performed the amplitude 
expansion for such a rotating cylinder of infinite radius 
and tT = 6.4 (corresponding to water), but, unfortunately, 
found that the coefficients remain close to the relax- 
ational line ci = — C3 for arbitrary values of the rotation 
rate B. Our results for the annulus, illustrated by the 



heavy line for (t = 6.7 in Fig. 1, confirm this for essentially 
all Prandtl numbers larger then 0.2. 
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FIG 1. Part of the phase diagram of the CGLE. Since 
the CGLE is invariant under ci — ci, C3 — C3, so is the 
phase-diagram; for convenience, we focus on the left half. 
The diagonal line ci — — C3 indicates the relaxational limit 
of the CGLE, while the curve C1C3 = 1 corresponds to the 
Newell criterion [Q]; outside of this hyperbola chaos occurs. 
The various spatio-temporal chaotic regimes and the lines 
Li, L2 and I/3 separating them are indicated, and the 
fat lines denote the results of our calculations. For a = 6.7, 
the coefficients stay close to the relaxational limit, similar to 
but for a around 0.15 the rotation-rate is able to tune 
ci and C3 from close to the relaxational limit to deep into 
the spatio-temporal chaotic regime. The variation with the 
dimensionless rotation rate Q is illustrated by the data-points 
for ^ = 1000 and J7 = 2000. 

In order to obtain an experimental realization of the 
CGLE with coefficients that can be tuned by the rotation 
rate over a much wider range, we have investigated con- 
vection in rotating annuli for a range of Prandtl numbers. 
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We take the circumference of the annulus large enough 
that the curvature can be neglected. This geometry dif- 
fers from a cylindrical one in that bulk-modes are quasi 
one-dimensional and can be described by the same am- 
plitude equations as the wall- modes, and that there are 
two wall-modes, localized near both side-walls. The in- 
finite radius limit allows the use of Euclidian instead of 
cylindrical coordinates and circumvents the wavenumber 
discretization that occurs in a finite container |^,^. In 
a system of finite radius, the two wall-modes will differ 
slightly, but we will ignore this. The main result of our 
calculations, illustrated in Fig. 1, is that for 0.1 ^ct^ 0.2 
traveling waves in a rotating annulus are predicted to be 
an experimental realization of a CGLE whose coefficients 
scan through the spatio-temporal chaotic part of the phase 
diagram when the rotation rate is changed. 

The relevant mode in the interesting low Prandtl num- 
ber regime is, however, not the analog of the wall-mode 
studied by Kuo and Cross |^ . We have found that the a- 
dependence of the wall mode is weak, and its coefficients 
remain close to the relaxational limit. However, in an an- 
nulus of finite width, the oscillatory bulk mode identified 
already long ago for an infinite container , becomes 

the primary bifurcation from the non-convective regime 
for a range of the rotation rates and 0.1 ^cr< 0.2. It is 
this bulk mode whose coefficients have the desired behav- 
ior. For even smaller values of cr, several nearby branches 
compete and no simple picture emerges |]l^ . 

Our calculations are based on standard methods to 
calculate the coefficients of the amplitude equations (see, 
e.g., Q), but are complicated by computational difficul- 
ties and the occurrence of several competing branches of 
solutions for a. Since a detailed account of the calcu- 
lations and the results is given in jlj], we confine our 
presentation here to a description of the basic setup and 
a summary of the most important predictions. 

The amplitude equations describing the slow modu- 
lations of a right-traveling mode with x, t-dependence 
^i{-kax+u)at) amplitude A, coupled to a left-traveling 
mode e^ii^ix+ujct) T^[i]^ amplitude B, are ||l[: 



Toidt + Vgdx)A 
To{dt 



vgd, 



e(l + icf))A + ^2(1 + ici)d^A - 
go{l - ic3)\A\^A - .92(1 - ic2)\B\^A , 

B = e{l + ico)B + e^(l + ici)dlB - 
go{l-iC3)\B\^B-g2il-iC2)\A\^B . 



(2a) 
(2b) 



For Raylcigh-Bcnard convection, e:—{R — Rc)/Rc, with 
R the Rayleigh number gaATd^ / kv, where g is the gravi- 
tational acceleration, a the thermal expansion coefficient, 
AT the temperature difference between bottom and top 
plate, d the height of the layer, n the thermal diffusivity 
and v the kinematic viscosity; Rc is the critical Rayleigh 
number. These quantities arise in the equations of mo- 
tion that describe the fluid-system in the co-rotating 
frame, which are the Navier-Stokcs equation with ad- 
ditional centrifugal and Coriolis force, supplemented by 
the heat equation and the mass conservation law |lQ] . As 



usual iP-pr]] , we apply the Boussinesq approximation and 
neglect the centrifugal forces. 

Below we shall summarize our findings for all the co- 
efficients and parameters in the amplitude equations. 
Since we find that g2 > go for the bulk-mode, the stand- 
ing waves are suppressed and the relevant dynamical 
states are traveling wave states with, e.g., Ay^O, B = 0. 
Upon rescaling time, space and the amplitude, Eq. ( |2a| ) 
then reduces to the CGLE (|^) in the frame moving with 
the group velocity Vg. 

As length-scale we choose d, so the top and bottom 
plates are at z = and z = l. We focus here on the case 
that the width of the channel is 1; for nearby values of the 
width similar behavior occurs, as detailed in |12|. The 
rotating (infinite radius) annulus is therefore character- 
ized by two parameters. The first is the Prandtl number 
(T v/k. The second is the dimensionless rotation rate 
VL:—Q,Dd^/v, where il^i is the angular velocity. In a typ- 
ical experiment, a is fixed, and the rotation rate f2 can be 
adjusted over a certain range up to values of order IC*. 
To be able to separate the hydrodynamic equations, we 
assume slip boundary conditions on the top and bottom 
plates, as in |^,|ll|. On the vertical side- walls we apply 
stick boundary conditions, which damp the mean-flow 
that plays a role for low Prandtl number convection with 
slip boundary conditions 

Vx = Vy = = dyO on y 0, 1 , (3a) 
dzVx = dzVy = Vz = = on z = 0,1 , (3b) 

where 9 is the deviation of the temperature from the 
conductive state proflle. It should be noted that our ver- 
sion of the stick boundary conditions is slightly simpler 
than those used in ||^. We have normalized the ampli- 
tudes such that \A\ represents the ratio of convected to 
conducted heat; the value of go therefore determines the 
so-called Nusselt number 

Bifurcation structure and linear stability. From the 
equations of motion we have determined i?c, which is the 
value of the Rayleigh number where convection sets in, 
and the corresponding critical wavenumber kc and fre- 
quency LOc as a function of a and fi. An important fea- 
ture of our system is that there exists, in particular for 
small fj's, a multitude of solutions to the linearized equa- 
tions, but only the mode with the lowest critical Rayleigh 
number is relevant. For a finite cylinder, these branches 
are discussed in detail in [|j . 

The relevant features of the bifurcation structure are 
illustrated in Fig. 2 and can be summarized as follows. 
(a) The linear stability analysis for the stationary modes 
does not dependent on the value of cr ||l^. (b) For all 
(7 > 0.2, the wall-mode is relevant for sufficiently large 
rotation rates, while a stationary (non-oscillatory) mode 
is relevant for small Q, — see Fig. 2(c). For cr = 6.7 the 
crossover between stationary and oscillatory mode occurs 
at « 27.5; the situation for (7 = 6.7 is representative for 
the whole range cr ^ 0.2. (c) The bulk-mode exists for 
values of cr that are comparable to where the Hopf bifur- 
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cation in an infinite layer occurs |I^jri|], and approaches 
this mode rapidly upon increasing the width of the an- 
nulus. (d) For 0.1 cr ^ 0.2 we find a band of rotation 
rates flmin <^ < ^max for which the bulk-mode is rele- 
vant. See Fig. 2(c), where f2mm = 140 and ^max = 5600 
are marked. The situation for cr = 0.15 studied below 
is representative, as Fig. 2(d) for a — 0.175 shows, (e) 
When a < 0.1, the situation becomes quite convoluted 
since additional modes become relevant for some range 
of Q's. 
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FIG 2. (a)-(c): The linear onset values of bulk- (BM) and 
wall-modes (WM) as a function of Q, for a — 6.7 and 0.15. 
In (c-d) we have rescaled the critical Rayleigh numbers by 
Q, to facilitate the comparison between the various modes. 
The stationary bifurcation is seen to be the primary bifurca- 
tion for small rotation rates. For all a > 0.2, the wall-mode 
constitutes the primary mode for sufficiently large Q as il- 
lustrated for the case a — 6.7. For a = 0.15, the bulk-mode 
exists for all fl > ^Imin ~ 140. Its critical Rayleigh number is 
smaller than the critical Rayleigh number of the wall-mode 
for ilmin < n < fimax ~ 5600. It IS in this range that the coeffi- 
cients of the CGLE can be tuned over a wide range. In (a) and 
(b) the corresponding critical frequencies and wavenumbers 
are plotted, (d): For values of the Prandtl number around 
0.15, the critical Rayleigh numbers of the wall-mode only 
show a weak dependence on the Prandtl number, whereas 
those for the bulk-mode strongly depend on a. As a result, 
the range of rotation-rates for which the bulk-mode is relevant 
strongly depends on a; for (t = 0.175, the values of flmin and 
^max are approximately 200 and 2500, while for a = 0.125 
(not shown) they are 122 and 10500. 

Amplitude expansion. The amount of computer time 
that is needed for the calculation of the nonlinear coeffi- 
cients go, C3, (72 and C2 of the coupled amplitude equations 
(^ is substantial; therefore, we can not scan all system 
parameters simultaneously. Instead, we have performed 
a "trial and error" search in the (SI, (T)-space. We shall 
not exhaust the reader with the data thus obtained but 



concentrate on the wall- and bulk-mode discussed above 
and to a few values of a that are representative for the 
various ranges of the Prandtl number. The dependence 
of the coefficients ci and C3 of (||) for these two modes is 
illustrated in Figs. 1 and 3. 

For the wall- modes, wc find that the precise value of 
a is quite irrelevant; the coefficients ci and C3 of the am- 
plitude equations are always near the ci = — C3 line, as 
the full line in Fig. 1 shows for a — 6.7. When is suffi- 
ciently large, the wall-modes at the opposite sides of the 
boundary are decoupled and we recover the result of Kuo 
and Cross Of course, we have performed extensive 
searches in parameter space to search for more interesting 
behavior of the coefficients ci and C3 , but for all Prandtl 
numbers larger than 0.2, the wall- mode is relevant, and 
the behavior of the coefficients of the amplitude equa- 
tions (H) is very much like that for a — 6.7. 

For all system parameters that we investigated, 172 re- 
mains smaller than go for the wall- modes, and there- 
fore the convection occurs in two counter-propagating 
traveling waves, that are localized near the opposite y- 
boundaries |l[]. The coefficients of the amplitude equa- 
tions for the wall-modes that are relevant for small a 
and > ^Imax 1 are similar to the coefficients that we find 
for the wall-mode for higher Prandtl numbers, in that 
ci f« — C3. In the other small-cr regime where the wall- 
mode is relevant, i.e., f2 < f2mm, there appears to be a 
tiny regime, close to flmin, where the coefficients ci and 
C3 might move away from the dissipative limit; however, 
the numerics are not decisive here. 

From the point of view of the amplitude equations, the 
wall-modes do not have many interesting features, and 
therefore we will focus now on the bulk-modes that are 
relevant for small Prandtl numbers and Qmin <^< ^max- 

The coefficients of the amplitude equations (^) for the 
bulk-mode are shown in Fig. 3 for a = 0.15. We find 
that the coefficients ci and C3 can be tuned over a wide 
range by the rotation rate. For the system parameters 
that we consider, 32 > go [see Fig. 3(d)] and the left- and 
right-traveling modes suppress each other. The convec- 
tion patterns thus consist of a juxtaposition of patches of 
left- and right-traveling waves, and after careful adjust- 
ments one may have the convection exclusively consisting 
of either a left- or a right-traveling wave, which warrants 
a description with a single CGLE. Note that for such a 
single wave the value of C2 is immaterial, since terms of 
the form jylpi? or jSpyl are zero. 

At other values of the Prandtl number a in the range 
0-1 ^cr^ 0.2, the main effect of a change in a is through 
a change in ^min and Qmax the coefficients of the ampli- 
tude equations depend of cr, but this dependence is rather 
weak, in the sense that Q still makes the coefficients vary 
over a wide range. This is shown in Fig. 1, where the 
path that the coefficients trace in the parameter space of 
the CGLE is shown for three values of a. 

Can the interesting Prandtl number range 0.1<cr^0.2 
be accessed experimentally? Compressed gases have 
cr > 0.7 and typical liquids have even larger cr's, while 
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liquid mercury and gallium have (t = 0.025 and (7 = 0.005. 
However, Rayleigh-Benard convection in superfluid '^He- 
'''He mixtures is known to behave to a very good ap- 
proximation as a convecting liquid with a Prandtl num- 
ber which can be tuned continuously between 0.02 and 
about 1 [Q. Convection in rotating cells and flow visu- 
alization have recently become possible for such mixtures 
Ipzl , and so this system provides a unique route to probe 
the regimes of phase-chaos and defect-chaos experimen- 
tally, and to compare to theoretical predictions 
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FIG 3. The coefficients of Eq. (2) for the bulk-mode for 
a — 0.15. The coefficients Ci and in particular C3 have a strong 
dependence of the rotation rate (c). For Q — > flmax, ci and 
C3 are close to the relaxational limit. When the rotation rate 
is decreased C3 changes sign, and at f2 « 1050, the Newell 
criterion is reached (ciC3 = l). A further decrease of the rota- 
tion rate pushes the coefficients deep into the spatio-temporal 
chaotic regime as Fig. 1 shows. 

In summary, we predict that Rayleigh-Benard convec- 
tion in a large rotating annulus is an attractive experi- 
mental realization of a supercritical CGLE with tunable 
coefficients for a number of reasons: (i) The onset of con- 
vection can occur either via a stationary or a Hopf bifur- 
cation; in the latter case the mode can either consist of 
a single traveling bulk-wave, or two counter-propagating 
wall-modes. The rotation rate can be adjusted to study 
the competition between these states, in analogy to the 
study of the co-dimension-2 points that occur in binary 
liquid convection Q. (ii) The quasi one-dimensional ge- 
ometry of the system warrants a description in terms 
of one-dimensional amplitude equations, for both the in- 
trinsic one-dimensional wall-modes and the intrinsic two- 
dimensional bulk-mode. (Hi) The onset of convection oc- 
curs via a forward bifurcation; for backward bifurcations. 



like those in binary liquids, amplitude equations can at 
most give a qualitative description of the patterns, (iv) 
The underlying basic equations for this system, i.e., the 
Navier-Stokes equations, are considerably simpler than 
the basic equations for convection in liquid crystals or 
binary liquids. For instance, in the latter system, it is 
hard to decide which aspects of the experimentally ob- 
served chaos can be described by the quintic CGLE, 
and which aspects are connected to physically relevant 
effects that are not captured in an amplitude description. 

We have benefitted from correspondence with R. E. 
Ecke, M. C. Cross and P. Kolodner. 
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